Linear and non linear response in the aging regime of the ID trap model 
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We investigate the behaviour of the response function in the one dimensional trap model using 
scaling arguments that we confirm by numerical simulations. We study the average position of the 
random walk at time t w + 1 given that a small bias h is applied at time t w ■ Several scaling regimes 
are found, depending on the relative values of t, t m and h. Comparison with the diffusive motion in 
the absence of bias allows us to show that the fluctuation dissipation relation is, in this case, valid 
even in the aging regime. 



PACS numbers: 75.10.Nr, 05.20.-y, 02.50.-r 

The one dimensional trap model (hereafter denoted as 
1DTM) has been the focus of renewed attention, both 
in the mathematical community [1,2], and also using a 
physicist approach [3,4]. This model was proposed in the 
70's to describe the properties of one dimensional disor- 
dered conductors [5,6]. Very recently, the direct relevance 
of this model for the dynamics of DNA 'bubbles' under 
torsion was emphasized [7]. Although the "annealed" 
version of the model, in which a new trapping time is 
randomly choosen at each step, is now well documented 
[8-10], the full analysis of the quenched model consid- 
ered here remains challenging. It was shown in [3,4] 
that, besides interesting dynamical localization proper- 
ties, the 1DTM exhibits different time scalings, depend- 
ing on which correlation function is considered. To be 
more specific, the probability of not moving between t w 
and t w + 1 - called H(t w + t,t w ) in [3] - and the proba- 
bility of occupying the same site at t w and t w +t - called 
C(t w +t, t w ) - exhibit different scalings: n scales as tjt v w 
with v < 1, whereas C behaves as t/t w . Then a natural 
question arises in this context: what is the relevant time 
scale that governs the response function of the particle to 
an external bias? This question is interesting in the con- 
text of the physical applications mentioned above, and 
also in the context of the out-of-equilibrium Fluctuation 
Dissipation Theorem that was much discussed recently 
[11]. 

Let us first briefly recall the definition of the 1DTM. 
Consider a one dimensional lattice, and associate to each 
site i a quenched random variable Ei > 0, the energy 
barrier, chosen from a distribution p(E). One follows the 
evolution of a particle driven by a thermal noise at tem- 
perature T on the lattice. The particle has to overcome 
the energy barrier Ei in order to leave site i and reach 
one neighbour. This naturally leads to an escape rate 
Wi of site i given by an Arrhenius law uii — Tq e~ Ei ' T . 
Note that the transition rates are a coarse grained de- 
scription of the underlying Langevin dynamics, which is 
not explicitely described in this model. The quantity 
Tq is a microscopic frequency scale, that will be set to 
unity in the following. Once that particle has escaped 
the trap, it chooses one of the two nearest neighbouring 



sites, with probability q- for the left one and q + = 1 — g_ 
for the right one. Two particular cases have already been 
studied in details in the literature, namely the unbiased 
case g+ = 1/2 [3,1,2], and the fully directed case q + = 1 
[12,6,13]. 

In order to compute a response function, one has 
first to define an external field to which the response 
will be associated. The simplest choice is to consider 
small deviations from the unbiased case, i.e. to intro- 
duce a small bias h, independent from the site, such that 
q± = (1 ± h)/2. Note that this is equivalent, as long as 
linear response is obeyed, to the introduction of a uni- 
form small force field F, which transforms the transition 
rates q± into q± = q± exp(±_Fa/2fcT). Here a is the lat- 
tice spacing, which is set to unity in the following; the 
correspondence is then h — F/2kT. Physically, this force 
is due to the external electric field in the case of conduc- 
tors [5], or to an asymmetry of the DNA composition 
along the chain in the model considered in [7]. 

As far as the response is concerned, the two natural 
choices are the average position of the walk after the 
bias is applied, or the average probability current. Inter- 
estingly, a formal relation exists between the former and 
the latter, namely that the current is the time derivative 
of the average position. Switching on a small bias h at 
time t w and measuring quantities at a subsequent time 
t w +t, we define for a given sample of the disorder the av- 
erage position Xh(t,t w ) and the total probability current 
Jh{t, t w ) as: 
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with 4> n -i.n = W^_ 1 ^ n P„-i - W r ^ n _ 1 P n being the 
local current. Taking the derivative of Xh(t,t w ) with re- 
spect to t, one has: 

= ^2 n ~^f = ^2 n i^n-l,n - (t>n,n+l] (3) 
n n 

Assuming periodic boundary conditions on a lattice of 



1 



size L, and letting eventually L — > oo, one can show that 
the sum reduces to ^2 n <j) n -i,n, finally leading to: 



dxfi 



(4) 



Note that in the 1DTM, transition rates W^ n±1 depend 
only on n, so that Jh(t,t w ) = for h = 0, leading to a 
value of x(t, t w ) — x{Q, 0) in the absence of bias, for any 
given sample. Averaging over the disorder, it is clear that 
Eq. (4) is also valid for the averaged quantities (x)h{t, t w ) 
and (J)h(t,t w ). Therefore in the following we shall focus 
only on the average position after a bias is applied. 

In this section, we shall give some simple scaling argu- 
ments in order to predict the behaviour of (x)h(t, t w ) as a 
function of the three variables t, t w and h. Interestingly, 
non trivial regimes appear due to the fact that the limits 
h — > and t, t w — > oo cannot be inverted. Note that the 
case t w = was studied in [6], where it was shown that 
a non trivial crossover line appears in the plane h,l/t. 
Let us recall briefly these scaling arguments, since this 
will be useful in the following. It is convenient to in- 
troduce the typical number N of steps of the walk after 
time t, and to express both (x)h and t as a function of 
N. It is clear that (x)h — Nh; now considering the typ- 
ical number Af s of sites visited by the walk, Af s can be 
approximately written as the sum of a drift contribution 
and of a diffusive one: 



(5) 



Consider first the case Nh 3> vJV, corresponding to 
Af 8 <~ Nh. Given that the trapping times n are dis- 
tributed according to p(r) = /i/t 1+m , the sum of M 
independent variables Tk behaves as Y^k=i Tk ~ M 1 ^. 
Since each site is visited of order N/Af s times, t can be 
expressed as: 



XT S ~ Nh + VN 



N ^ 

•=-jV,/2 



NM S " (6) 



which can be rewritten as N 
reads: 

(x) h ~ h"t* 



t^h^- 1 . Finally, (x) h 



(7) 



The criterion Nh 3> translates into h^t^ ^> 

ft ( M -i)/2 tM /2 j or equiva iently t > t h ~ fc-a+JO/M. On 

the contrary, if Nh < VN then Af s ~ \/N, and N and t 
are related through N <~ t 2fiu , so that: 



(x) h ^ht 2 ^ 



(8) 



with v = 1/(1 + fi). Therefore, the response is linear in 
h but non linear in t for t -C th, and non linear both in 
h and t for t » th- 

Let us now turn to the response in the aging regime 
t w » 1. Then if t is small enough (to be specified later), 



the walk will essentially evolve within the space region 
of size £(t w ) already visited (the diffusion correlation 
length). The average trapping time r(t w ) within this 
region is defined by: 



r(t w ) 



j; 



fidr 
T i+n 



t 



(9) 



which is the average value of t computed from the dis- 
tribution p(r), taking into account the natural cut-off in- 
duced by the dynamics t w [3]. Therefore, as long as the 
particle does not escape from the initial region, the aver- 
age position should drift at a constant velocity h/r(t w ) 
(the lattice spacing a is taken as the unit of length), lead- 
ing to: 



(x) h ~ ht/t 



"(l-M) 



(10) 



This short time regime is limited by two conditions: first, 
\€(t w + 1) — i{t w )\ <C £(*«,), which implies t <C t w , and 
also \(x)h \ <C £,(t w ), which requires (<f, where t* is a 
new time scale defined by t* = t w /h. Note however that 
this time scale is only relevant if t* <t w . If we are in the 
opposite limit t 3> t w , then t w no longer plays any role 
and one recovers the results found above in the particular 
case t w = 0. So one has to distinguish between several 
regimes, depending on the relative values of t, t w and t*: 

• t* » t w (or h <C h* <~ tw'"'")- I n tn i s case : three 
different regimes appear: 



(x)h ~ ht/t 
(x) h ~ht 2 » v 
(x) h ~h»t» 



f(i-f*) 



t < i w 
i w < t < t h 
* > th 



(11) 
(12) 

(13) 



where th ^ h as above. 



• i* <C tu, (or /i > h*). The behaviour of (a;)/, can 
be summarized as follows: 



(x)h 



h" f 



(14) 
(15) 



It is interesting to reformulate the above equations in 
terms of scaling functions: 



{x) h = hty h{t/t w ) h^h* 
(z)h=C Hht/t w ) h^h* 



(16) 
(17) 



with the following asymptotic behaviour for the functions 
fi(z) and f 2 (z): 

h(z)~z fi(z)~z 2 »" (z»l) (18) 

f 2 (z)~z (*«1), / 2 (*)~*" (z»l) (19) 
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The scaling function fi(z) only accounts for the time 
regimes described by Eqs. (11,12), i.e. for times smaller 
than th- However, since we are in the small bias case 
h <C t~^ v , th is much larger than t w , so that this 
crossover scale is difficult to evidence numerically. 




t/t 

w 

FIG. 1. Rescaled response in the small Geld regime 
(h < h*(t w )), using Eq. (16), for h = 2.KT 3 , 5.KT 3 , 8.1(T 3 
and t w = 10 3 , 10 4 , 10 5 (p = 1/2). The resulting collapse is 
good, even though some deviations appear both at short and 
large times. The predicted asymptotic scaling behaviour is 
well obeyed (dashed lines). 

We have not found the way to include all the t, t w and 
h regimes in a single scaling function. We now report 
on the numerical results on the response function, and 
compare them to the scaling predictions of the previous 
section. In the small field regime h -C h*(t w ), the scaling 
predicted by Eq. (16) is well satisfied. Fig. 1 shows the 
resulting collapse of the data for three different values of 
h (h = 2.1CT 3 , 5.1CT 3 and 8.1CT 3 ) and three different 
values of the waiting time t w (t w = 10 3 , 10 4 and 10 5 ) at 
temperature [i=\- 

In the opposite regime, h » h*(t w ), Eq. (17) is also 
well obeyed, as shown on Fig. 2 for bias h = 0.1, 0.2, 
0.3, waiting times t w = 10 4 , 10 5 , 10 6 , 10 7 and fj, = \. 
Note that for clarity, data corresponding to h = 0.2 and 
0.3 are presented for t w = 10 5 only. We have checked 
for several values of /j, that the short time and large time 
behaviour of the scaling functions fi(z) and fziz) given 
by Eqs. (18,19) are also correctly predicted - see Fig. 1 
and Fig. 2 for the case /i = §• 

If one wishes to draw a link with the correlation func- 
tions C(t w +t,t w ) and U(t w +t,t w ) defined above (see 
also [3]), it might appear that (x)h(t, t w ) should be asso- 
ciated to C(t w +t, t w ), due to the t/t w scaling in Eq. (16). 
However, this is limited to the very small bias case. Inter- 
estingly, in the opposite regime h h*(t w ), (x)h(t, t w ) is 
a function of t/t* (t* ~ t^/h), up to a prefactor depen- 
dent on t w . So at fixed h, (x)h(t,t w ) scales as t/t v w , as 
U(t w +t, t w ) does. In fact, this relation is not purely for- 



mal, but indeed corresponds to the underlying physics. 
During the unbiased time interval [0, t w ], the walk typi- 
cally visits deep traps with characteristic time ~ t v w . 




FIG. 2. Rescaled response in the large bias regime 
(h > h*(t w )), using Eq. (17), for h = 0.1, 0.2, 0.3, t w = 10 4 , 
10 5 , 10 6 , 10 7 (see text) and /i = 1/2, showing a very good col- 
lapse of the data, except for small t where finite time effects 
become noticeable. The asymptotic behaviour of the scaling 
function is also well predicted (dashed lines). 

Once the bias is applied, the evolution is dominated 
by these deepest traps, before eventually reaching a long 
time regime independent from t w . If the bias is very 
small, the walk will visit these traps a large number of 
times (of order t^f [3] ) , and the aging dynamics thus re- 
sembles closely that of C(t w + t,t w ) in the absence of 
bias. On the contrary, if h is large enough, the walk will 
visit only a few times (of order 1/h) the deepest traps 
occupied at time t w , so that the aging dynamics is dom- 
inated by the time needed to leave the deep traps for 
the first time, in close analogy with U(t w + t,t w ). More- 
over, it is worth mentioning that if one applies the bias 
from t = instead of t w , and compute the resulting aver- 
age displacement between times t w and t w + t, different 
scalings are obtained. For small fields, the displacement 
(x)(t,t w ) behaves as (x)(t,t w ) = ht 2 Jf v fi(t/t w ) as in the 
previous case, but for larger h, one finds a scaling of the 
form (x)(t, t w ) = t^f 3 (ht/t w ). So in this case, a t/t w 
scaling also appears, but for a different reason: because 
of the bias, the trapping times reached after a time t w are 
now of order t w itself instead of t v w , but the walk visits 
these traps a finite number of times (~ 1/h) after time 
t w . 

Besides the correlation C and n considered hereabove, 
the natural self-correlation associated to (x)h is the 
mean square displacement restricted to the time inter- 
val [t w ,t w +t], in the absence of bias field: 

(Ax 2 )o(i, t w ) = ([x(t w +t)- x(t w )] 2 ) (20) 

where the brackets (• • •) denote both a thermal av- 
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erage and an average over the disorder. Using the 
effective trapping time f(t w ), one can also estimate 
(Ax 2 ) (t,t w ) in the short time, diffusive regime as 
{Ax 2 )o(t,t w ) <~ t/r(t w ). This expression can only be 
valid if (Ax 2 ) (t, t w ) <C ^(t w ) 2 , which yields t <C t w . In 
the opposite regime t ^> t w , one recovers the t w = 
result (Ax 2 ) {t,t w ) <~ £(£„ + t) 2 — £,(t) 2 . In summary: 

(Az 2 )o(M™) ~ i/C^ *«*,» (21) 
(Ax 2 ) (t,t w ) ~t 2 ^ t^>t w (22) 

This means that (Ax 2 }o(t, t w ) can be written in the scal- 
ing form: 

(Ax 2 )o(t,t w )=t 2 w ^g(t/t w ) (23) 

Fig. 3 displays the numerical data obtained for 
(Ax 2 )o(t, t w ), for waiting times t w ranging from 10 3 to 
10 7 , using Eq. (23). One can see that the collapse is 
perfect. 
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FIG. 3. Rescaled diffusion {Ax }o(t,t w ) in the absence of 
bias, using Eq. (23), for t w = 10 3 , 10 4 , 10 5 , 10 6 and 10 7 , and 
fi — 1 /2. The different curves are indistinguishable to the eye. 
The dashed lines indicates the slopes 1 and 2/3 respectively. 
Note however that for computational reasons, time t is lim- 
ited to 10 7 . Inset: FDR plot {x)h(t,t w )/h versus (Ax 2 )o(t,t w ) 
fort w = 10 3 , h = 5.1(T 3 andu = 1/2. FDR is clearly satisfied 
over the whole range of time. 

Given the quantities computed above, it is natural to 
test the Fluctuation-Dissipation (or Einstein) relation 
(fdr). This is done in the inset of Fig. 3, which dis- 
plays (x)h(t,t w )/h versus (Ax 2 }o(t, t w ) in log-log scale, 
for h = 5.1CT 3 , t w = 10 3 and /z = 1/2. This relation 
appears to be very well satisfied over the whole range of 
time, although the system is strongly out of equilibrium. 

We now give a general argument in order to demon- 
strate the validity of the FDR for the trap model in the 
aging regime. It was shown in [6] that for a given config- 
uration of the disorder and a given initial position x(t w ), 
the following fluctuation dissipation relation holds, in the 
limit h — ► 0: 



(Ax) h = (Ax) + h[(Ax 2 ) - (Ax) 2 ] (24) 

with Ax = x(t w + t) — x(t w ), and t is finite. For the pur- 
pose of clarity, we distinguish here between average on 
thermal histories after t w ((•••)) and average over ther- 
mal histories before t w - thus over x(t w ) - and over the 
disorder (~) . Now applying this second type of average 
to the previous equation, we get: 



(x) h (t,t w ) = —{Ax 2 ) (t,t w ) (25) 

where we have used the fact that for the trap model 
{Ax)o = - see Eq. (4) -, and h has been replaced 
by its 'physical' expression F/2kT. So we conclude that 
the fluctuation dissipation relation is indeed valid in this 
out-of-equilibrium and disordered system, with a tem- 
perature equal to the bath temperature T. In particular, 
this implies that the scaling functions /i(.) and g{.) arc 
identical. Note that such an "aging Einstein relation" 
has already been found in the annealed version of the 
model [9]. 

As a conclusion, we note that the influence of an exter- 
nal bias in disordered system can be highly non trivial. 
In the simple one dimensional trap model discussed here, 
we already find several regimes in the t,t w ,h 'cube'. We 
have extended the above arguments to the Sinai model 
with an external bias (for a review, see [14]). Already for 
t w — and small external force F, one finds in general 
four different time regimes for the average displacement 
(x). There is in particular a regime where (x) grows as 
In 2 t, but with an F independent prefactor, before the 
asymptotic regime where (a;) <~ t aF sets in. In other ex- 
amples, such as a walker on a percolation network, the 
response can even be non monotonous with F [6,15]. 
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